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Abstract 

QCD is expected to be in the color-flavor locking phase in high baryon density, which 
exhibits color superconductivity. The most fundamental topological objects in the color 
superconductor are non-Abelian vortices which are topologically stable color magnetic flux 
tubes. We present numerical solutions of the color magnetic flux tube for diverse choices 
of the coupling constants based on the Ginzburg-Landau Lagrangian. We also analytically 
study its asymptotic profiles and find that they are different from the case of usual super- 
conductors. We propose the width of color magnetic fluxes and find that it is larger than 
naive expectation of the Compton wave length of the massive gluon when the gluon mass 
is larger than the scalar mass. 
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1 Introduction 



One of the important problems for understanding the strong interaction is to determine the 
phase diagram of QCD. In the very high density region with a large chemical potential at low 
temperature, it is expected that QCD is in the color-flavor locking (CFL) phase which exhibits 
color superconductivity [TJ [2]. There is the S77(3)l x S77(3)r flavor symmetry acting on light 
quarks u, d and s when their masses are very small compared to the chemical potential. The 
S77(3)c color symmetry is locked with the flavor symmetry by the condensations 



where ^l.r are the quark fields. Apart from the discrete symmetry the symmetry of the system 
is 



where the first element and the rests stand for gauge and flavor symmetries, respectively. We 
assume that the ground state is the positive parity state $r = <3>l = $ which would be determined 
by the instanton effect, so the chiral symmetry is broken to S"?7(3)l+r which we write SU(3)p. 
The color superconductivity is expected to be realized in the core of a neutron star. 

In the usual superconductors the gauge group £7(1) for the electro- magnetic force is broken in 
the vacuum where electrons condensate to make a Cooper pair. When the superconductor is in 
the external magnetic fields which are above the critical value, magnetic fields inside the supercon- 
ductor are squeezed and quantized. Consequently there appear Abrikov-Nielsen-Olesen(ANO) 
vortices as magnetic flux tubes which topologically wind the spontaneously broken 17(1) OH]. 
The stability of the superconductor under the external magnetic fields is determined by static 
forces among vortices. If the force is repulsive (type II) the superconductor is stable whereas if 
the force is attractive (type I) it is unstable. In the former the vortices constitute the so-called 
Abrikosov lattice [3]. Therefore understanding the interaction between vortices is very important 
to study stability of superconductors. 

The types of color superconductor (type I or type II), which characterize a response to external 
magnetic fields, were studied in [3] by considering a domain wall separating the normal phase 
and the superconducting phase. They concluded that it is of type I in the weak coupling region. 



(1.1) 



G = SU(3) C x SU(3) L x SU(3) R x U(l) 



(1.2) 
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It is, however, more important to study vortices in color superconductor, in order to understand 
its stability or a response to (external) color magnetic fields and/or rotation of the media which 

><fl 



exhibits superfluidity toco Color magnetic vortices were studied in (6J [3 [8] but these are unstable 
to decay because they are not topologically protected. On the other hand the superfluid vortices 
which wind the broken £/(1)b appear in a response to rotation 0, U\- They are topologically 
stable but are dynamically unstable as we discuss below. 



The most fundamental objects are non-Abelian vortices carrying color magnetic fluxes [TOl 
[TTT IT2] . which were called the semi-superfluid vortices. Let us explain these objects. The order 
parameter space, parameterized by (would-be) Nambu-Goldstone(NG) modes of the symmetry 
breaking in the CFL phase, is 

^3 

These NG modes are eaten by the SU(3) gauge fields (gluons) except for the U(l) part. This 
has a non-trivial first homotopy group; 

TTi (M) = TTi [U(3)} = Z. (1.4) 

The U(1)b vortex found in [HI E] is of the form $ = v /(r)l 3 with the boundary conditions 
f(r—> oo) = 1 and f(r — 0) = 0. However this does not have minimum winding number but has 
the triple of the minimum element. The fundamental vortex-string with the minimum number 
is a non-Abelian vortex [TUl dH d2] , for instance given as 

$ = v diag(/(r)e ie , g(r),g(r)), ~ ^0- [1 - h(r)} diag(2, -1,-1) (1.5) 

with the boundary conditions f{r^ oo), g(r oo) = 1, h(r — > oo) = and f(r — 0) = 
0, g>'(r = 0) = 0, /i(r = 0) = 1. The asymptotic behavior (r — > oo) of the scalar field of the 
non-Abelian vortex is 

$ -> t;diag(e ie , 1, 1) = ue^diag (e*S e , e _i i e -<! ^ ^ 1 3 , (1.6) 



1 In a realistic situation such as a neutron star, the external color magnetic field cannot be considered because 
it is in the confining phase outside the core of a neutron star. Our statement about a response to color magnetic 
field is of a theoretical interest. The vortex here should be created under a rapid rotation of superfluids, or under 
a phase transition in a rapid cooling of a neutron star by the Kibble-Zurek mechanism [10] . 
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where the last denotes a gauge transformation by U(r,6) = diag l l eE ( r ) ^ e z l E ( r ) ^ e l l E ( r )^j vvit 



ifl 



an arbitrary function E(r) satisfying the boundary conditions E(r = 0) = and E(r — > oo] 
From Eq. (11.61) one can understand that the non-Abelian vortex has the U(1)b winding 1/3 of 
that of the U(1)b superfluid vortex. The U(1)b symmetry is global and so the tension (the energy 
per unit length) of a vortex-string is logarithmically divergent in infinite space. Since the tension 
is proportional to the square of the winding number, the non-Abelian vortex has the tension 1/9 
of the one of the U(1)b superfluid vortex. One U(1)b superfluid vortex has the triple amount of 
the sum of the energies of three well-separated non-Abelian semi-superfluid vortices. Therefore 
the decay of U(1)b superfluid vortex into three non-Abelian semi-superfluid vortices is inevitable 
from a point of view of the energetics. 

The non-Abelian semi-superfluid vortex (jl.5p carries a color magnetic fluxjfl Color magnetic 
fluxes also exist in quark gluon plasma [H] but they are unstablejj Contrary to those, color 
magnetic flux tube (jl.5p is topologically (and dynamically) stable. The vortex solution (11.51) 
breaks the color-flavor locked symmetry S77(3)c+f down to its subgroup [SU(2) x {7(1)] C+F . 
Consequently there appear further NG zero modes [11] 

CP 2 = ^ (3)c+F (1 7) 

[SU(2)xU(l)] c+F 1 • ] 

This space parameterizes a continuous family of the vortex solutions. These zero modes are 
called orientational zero modes|] All solutions of the continuous family in Eq. (11.71) have the 
same tension and the same boundary condition (II. 6p up to a regular gauge transformation. 
Therefore the orientational zero modes of CP 2 (I1.7P are normalizable and can be regarded as the 
genuine moduli (collective coordinates) of the vortex [TB|. It corresponds one-to-one to the color 



2 This transformation is well-defined in the whole space because of the triviality of the first homotopy group: 
m[SU(N)] = 0. 

3 The authors in [13] studied different kind of magnetic flux tubes in dense QCD. Although the authors called 
them "gluon vortices" , those are magnetic for usual electro-magnetic force of [7(1)em (mixed with one color 
component of SU(S)c)- The important difference between their flux tubes and our color magnetic flux tubes is 
that their vortices arc generated by the unbroken generator and therefore topologically trivial, while our vortices 
are generated by the broken generators and consequently topologically non-trivial. We will give some comments 
in the discussion. 

4 Those flux tubes are topological unstable but it was discussed in [15] that they are dynamically metastable. 

5 The idea was brought from non-Abelian vortices in supersymmctric U (N) gauge theories. In this case the 

overall phase E/(1)b is also gauged and the vortices are local vortices with finite tension [16], unlike non-Abelian 
semi-superfluid vortices. See [17] for a review. 



3 



magnetic flux which the vortex carries. 

As in usual superconductors, the (in) stability of the color superconductor in the presence of 
the (external) color flux is determined by the interaction between non-Abelian semi-superfluid 
vortices (color magnetic fluxes) given in Eq. (11.51) . The asymptotic interaction between two well- 
separated non-Abelian semi-superfluid vortices has been calculated pj], [12] in which the universal 
repulsion has been found. This calculation is valid when the distance between them is much larger 
than the Compton wave lengths of massive particles, which are essentially the penetration depth 
and the coherence length. In this region, semi-superfluid vortices are essentially U(1)b global 
vortices with the winding number 1/3 as in Eq. ( 11.61) and, in fact, the static force between them is 
1/3 of that between two U(1)b vortices [TTl[T2] . This result implies important consequences. First 
the color superconductor is stable in the presence of the (external) color magnetic fields. Non- 
Abelian semi-superfluid vortices will constitute a lattice, at least when the lattice spacing is much 
larger than the penetration depth and the coherence length. Second the U(1)b superfluid vortex 
is unstable to decay into three non-Abelian semi-superfluid vortices because of the repulsion 
among them, at least for large fluctuations. Each semi-superfluid vortex carries different color 
magnetic flux with the total color being cancelled out. In the core of neutron star, one can expect 
that the U(1)b superfluid vortices are first created in a response to rotation of the star. Then 
each of them must be broken into three semi-superfluid vortices which will constitute the lattice 
of color magnetic flux tubes. Phase boundary of the CFL and hadronic phase was studied in [19] 
in an application to the neutron star physics. 

However the interaction obtained in [TTl [T2"] is not valid when two vortices are closer such 
that the distance between them is of the order of the Compton wave lengths. We need to know 
the core structure of the vortex in order to study short range interactions. In the analysis of 
Balachandran et. al [TU], they used an approximation for profile functions in (jl.5|) : they assumed 
a constant g and solved equation for / only. This is not a good approximation around the core 
of the vortex. In fact, as we show in this paper, the profile of g significantly differs from the 
constant for some parameter region. 

In this paper we study the profile functions of non-Abelian semi-superfluid vortices in detail 
as a first step to study the short range interaction of them. First we study the asymptotic tails 
of the profile functions by analytically studying the equations of motion, and find that they 
are different from the case of the ANO vortex in usual superconductors where the tails of the 
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U(l) gauge field with mass m e and the scalar field with mass rnn decay exponentially as e~ me 
and e~ mH , respectively; on the other hand the both tails of the scalar and gauge fields of a 
non-Abelian superfluid vortex behave as e~ m with the lighter mass m among the masses ma for 
massive gluons and m x for massive traceless scalar fields. We then construct numerical solutions 
for diverse choices of the coupling constants, by using the relaxation method with the appropriate 
boundary conditions. By the numerical solutions we determine the coefficients in the asymptotic 
solutions. We propose the width of the color magnetic flux by the diffract iveness weighted by 
the magnetic flux. We calculate it numerically and confirm the above estimation of the tails. 
Therefore we conclude that the width of the color magnetic flux is significantly different from 
the naive expectation of the Compton wave length (the penetration depth) of the massive gluons 
when the gluon mass mg is larger than the mass m x of the traceless scalar fields; In general the 
width approaches to a constant depending of m x not behaving as m^ 1 . 

Throughout this paper we turn off the electro-magnetic interaction C/(1)em i n order to study 
purely non-Abelian aspects of vortices. The inclusion of it does not change the properties found 
in this paper. However the electro-magnetic interaction explicitly breaks the flavor symmetry 
S77(3)f, since {7(1)em is embedded into S77(3)f- Consequently the orientation modes of CP 2 get 
masses. This aspect remains as a future problem. 

This paper is organized as follows. In Sec. [2] the Ginzburg-Landau Lagrangian is studied. We 
explain the symmetry structure in detail especially paying attention to the discrete symmetries. 
We also give the mass spectra of the CFL vacuum. In Sec. [3] we study the equations of motion 
in the cylindrical coordinates in the presence of a vortex. In Sec. H] asymptotic behaviors of the 
profile functions are studied analytically. In Sec. [5] we give numerical solutions of the profile 
functions and determine the coefficients of the asymptotic profile functions. We numerically 
evaluate the width of the color magnetic flux and compare it with the Compton wave length 
of the massive gluons. Sec. [6] is devoted to conclusion and discussions. We discuss the force 
between two semi-superfluid vortices at the short distance. We discuss that the inter-vortex 
force mediated by the exchange of the massive particles becomes comparable with the inter- 
vortex force mediated by the exchange of the massless U(1)b NG boson. In Appendix we discuss 
general ansatz of the vortex in the diagonal entries. 
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2 The Non-Abelian Landau- Ginzburg Model 



Our starting point is the Ginzburg-Landau Lagrangian [20 
C = Tr 



where the last constant term is introduced for the vacuum energy to vanish. Our notation is 
Vi = di - ig s Ai, Fa = diA 3 - djA, - ig s [A h A,} and Tr[T a T b ] = 5 ab for a = 1, 2, • • • , 8. g s is 
S77(3)c gauge coupling constant. For the stability of vacua, we consider the parameter region 
yU 2 > 0, A2 > and 3Ai + A2 > 0. The action of color, flavor and baryon symmetries on $ is 
given by 

$ -> e ie U c ®U F , U c G SU(3) C , U F G SU(3) F , e w G U(1) B . (2.2) 

There is some redundancy of the action of these symmetries. The actual symmetry is given by 

c _ SU(3) c xSU(3) F xU(l) B ^ 

(Z 3 )c+B X (Z 3 ) F+ B 

where the discrete groups in the denominator defined as follows do not change <3> and are removed 
from G; 

(Z 3 ) c+ B : (w k l 3 , 1 3 , w~ k ) G SU(3) C x SU(3) F x U(1) B , (2.4) 
(Z 3 ) F+B : (1 3 , w k l 3 , w~ k ) G SU(3) C x SU(3) F x U(1) B , (2.5) 
w = e 2m/3 , fc = 0,l,2. (2.6) 

For later use let us redefine the discrete symmetry in the denominator as (Z 3 ) c+B x (Z 3 ) F+B — 
(Z 3 ) c+F x (Z 3 )c_ F+B with 

(Z 3 ) c+F : (w k l Nl w~ k l N , 1) G SU(3) C x SU(3) F x U(1) B , (2.7) 
(Z 3 )c„ F+B : (w k l N , w k l N , w~ 2k ) G SU(3) C x 5C/(3) F x [7(1) B . (2.8) 

Next we discuss symmetry breaking in the vacua. By using the symmetry G, one can choose 
a vacuum expectation value (VEV) as 

<*> = Bl * 1,2 s ^atW > (2 ' 9) 



The stiffness parameter in front of the kinetic term of <J> is set to be 1 by a field redefinition. 



without loss of generality. By this condensation the gauge symmetry SU (3)c is completely 
broken, and the full symmetry G is spontaneously broken down to 

R _ SU(3) c+ f x (Z 3 )c-f + b _ SU(3) C+F (21Q) 

(^3)c+B X (^3)f+B (Z^C+F 

with 

SU(3) C+F : (U, U\ 1) G SU(3) C x SU(3) F x [7(1) B . (2.11) 
Therefore as denoted in Eq. (11. 3p the order parameter space (the vacuum manifold) is given by 

M=G/H= smc^m.^ m (212) 

(^3)c-F+B 

with 

S77(3) C -f ■ (U, U, 1) G SU(3) C x 5C/(3) F x [7(1) b . (2.13) 

This space is parameterized by SU(3) would-be NG bosons, which are eaten by S77(3)c gauge 
bosons (gluons), and one massless NG boson of the spontaneously broken U(1)b- 

The mass spectra around the Higgs vacuum (12. 9p can be found by perturbing $ as 

$ = ,l 3 + ^l3 + ^V. (2.14) 

The trace parts <p and ip belong to the singlet of the color-flavor locked symmetry whereas the 
traceless part x an d C a belong to the adjoint representation of it. The SU (3) gauge fields (gluons) 
get mass with eating £ a by the Higgs mechanism. The masses of fields are given by 

m 2 G = 2 g y, ml = 2fx 2 , mj = 0, m\ = AX 2 v 2 , (2.15) 

where mc is the mass of the SU(3) massive gauge bosons (gluons) and ip is the NG boson 
associated with the spontaneously broken E/(1)b symmetry. The trace part <fi and the traceless 
part x °f ^ are massive bosons. At the weak coupling regime [7] these masses are shown to be 

m G 3> = 2m x . (2.16) 

We would like to stress that this system has four different mass scales. As in usual super- 
conductors, itIq 1 is the penetration depth and m~ l is the coherence length. On the other hand, 



7 We thank N. Yamamoto for this point. 



the massive boson with mass and the massless NG mode ip for the spontaneously broken 
£/(1)b are typical for the superfluid system with a superffuid vortex of the size m^ 1 . Therefore 
the system is mixed up with the (non-Abelian) superconductor and the superfluid. This is why 
the authors of pU] called this system the semi-superfluid. 



The interactions between two semi-superfluid vortices at large distance r 3> maxjm" 1 , m^ 1 , ttlq 1 } 
are interpolated by the U(1)b NG mode. As a result they repel each other by the long range 
force [TT] . 

Once the vortices are placed at the distance of the scale (9(m~^ G ), we cannot ignore exchange 

of the massive particles, and the interaction must depend on mass ratios 

m G _m x 
7i = , 72 = — ■ (2.17) 

m x 

Therefore, the type of relatively short range interactions would not be so simple, unlike the 
typel/II classification for usual superconductors. For instance, there is no critical coupling limit 
where the interactions are completely cancelled out in the case of color superconductor. The 
actual interactions at intermediate distance may be quite complicated. See Sees. @] and [6] for 
further discussion. 

For later convenience, let us rewrite the Landau- Ginzburg potential in terms of the dimen- 
sionful parameters fi 2 = mi/2, Ai = (mi — m x )/(12v 2 ) and A2 = m^/(4t> 2 ): 



v=^(^[* , *-" 2 i3]) 2 +^[<* , «> 2 ; 

where (A) stands for the traceless part of an N by N matrix A: (A) = A — (l/N)TrA. 
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3 Equations of Motion for a Vortex 



We would like to consider topologically stable vortex configurations supported by the first ho- 
motopy group in Eq. (jl.4p . Here we are interested in the minimally wound vortex solution. We 
make a diagonal ansatz for the minimally winding vortex in the cylindrical coordinates (r, 9, £3), 

/ e ie f[r) n ^ ' —9/3 n n \ 

g(r) 



$(Y, 



V 














9(r) ) 



1 ^ [1 - h{r)] 



9s 



\ 



-2/3 
1/3 
1/3 J 



(3-1) 



All other solutions are generated by the color-flavor symmetry. For later convenience let us 
rewrite the solutions (13.11) as 

*{r,B) = eKv^VTO (^) To _ y| G(r)T8 \ (3 . 2 ) 

1 enx> h 



g s r 2 L ' v /J V 3' 



Mr,0) = --^r\l-h{r)] x r-T, (3.3) 
with the U (3) generators 



T = i=diag(l,l,l), T 8 = -Ldiag(-2, 1, 1) (3.4) 

and the redefined profile functions 

F = f + 2g, G = f-g. (3.5) 

Note that the first term proportional to T in $ in Eq. (I3.2p is invariant under the color-flavor 
locking symmetry H in Eq. fl 2 . 1 [) while the second term proportional to T 8 (traceless part) breaks 
H down to U(2). This symmetry breaking by T§ gives rise to the NG zero modes SU(3)/U(2) ~ 
CP 2 associated with the vortices as in Eq. ( j 1.7ft pj]. Therefore the non-zero profile function 
G(r) plays a role of the order parameter for this breaking of the color-flavor symmetry. With the 
ansatz (13.31) . the color- magnetic flux is given by 

VEh' 



12 



3g s r 



T 8 . (3.6) 



Equations of motion for the profile function f(r),g(r), h(r) are of the form 



f + 7- i2h ^ l f ~ % (f 2 + ^ " 3 > " T f < /2 " & = (37) 

»" + 7 " ^» " (f + V - 3) + ^ 9 (f - 9 >) = 0, (3.8) 
h „ _ W_ _ ™| (j2(A _ 1} + /2(2ft + 1}) = Q (3 g) 

We solve these differential equations with the following boundary conditions 

(f,9, h) -»• (1, 1, 0) as r -> oo, 
(f,g',h)^ (0,0,1) as r^O. 

The third terms in the left hand side of Eqs. (13.71) and (13.81) are typical for global vortex config- 
urations which leads logarithmic divergence of the tension. 
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(3.10) 



The tension of the vortex-string is given by 



T 




(h-1) 



9r 2 



2 




+ v 



(3.11) 



V 



i, 2 [(/ 2 + 2, 2 -3) 2 m 2 +2( / 2 -, 2 ) 2 m 2 



(3.12) 



This diverges as T 



y-logL/r with L being an IR cutoff scale, the size of the system, and 



ro being a typical scale ro ~ . The factor 1/3 reflects the fact that the minimum winding 
vortex winds 1/3 of £7(1)b- 

We will see that the color-magnetic flux in Eq. ( 13.61) is well squeezed and becomes well localized 
tube even though the energy of the vortex itself logarithmically diverges. 

Here let us make comments on related vortices in the CFL phase. Another color magnetic 
vortex in the CFL phase studied in [8] is generated by only T 8 (mixed with electro magnetic 
U(1)em) without the part of T in Eq. (13. 2p . For the single- valuedness of fields the minimum 
winding is the triple of our color magnetic flux, and consequently it carries the triple amount 
of ours. It is, however, unstable because tti[SU(3)c} is trivial. On the other hand, the global 
C/(1)b vortex which is made by only the generator of global symmetry T has been studied in 
(9j[7j. It carries neither a color magnetic flux nor internal orientations. Since such configuration 
cannot be combined with SU(3), it is necessary to wind 2tt along To, and consequently it has 
the triple winding of our color magnetic flux. As stated in the introduction, it decays into three 
semi-superfluid vortices which have different color fluxes with the total flux cancelled out. 

Before closing this section, let us give a comment on a relation to a non-Abelian global vortex 
which appears when the chiral symmetry is spontaneously broken in the U(N) linear sigma 
model [22} [23} EH [25]. The linear sigma model is realized from our Lagrangian by just ungauging 
SU(N)c, namely turning off the gauge coupling g s = (rric = 0). At the level of equations 
of motion, it is enough to set h(r) = 1. Then Eq. (I3.9P becomes trivial and we are left with 
Eqs. (13.71) and (13. 8p (with h = 1) which are nothing but equations of motion for the non-Abelian 
global vortices [23} [25]. However h(r) = 1 is inconsistent with the boundary condition (13.101) . and 
so the non-Abelian global vortices cannot be obtained in a continuous limit of the non-Abelian 
semi-superfluid vortex. 
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4 Asymptotics 



Let us next study the asymptotics in the region r ^> max{m~ 1 , m^ 1 , m^ 1 } . The configuration is 
almost vacuum in such region, so it is useful for us to use the profile function F, G rather than 
f,g in Eqs. ( 13. 7ft and ( 13.81) . Let us perturb the fields by 

F = 3 + 5F, G = 5G, h = 8h, (\8F\, \8G\, \8h\) < 1, (4.1) 

and linearize the equations of motion as 



A -'4-9^)*G= Jj», (4.3) 

Sh" - — - mlSh = -m 2 c SG, (4.4) 
v 3 



where A = ^d r (rd r )u The last equation can be rewritten with 5h = Sh/(mG^) as 

1 \ 2 6G 
A-m 2 G - — ) 5h = -m G — . (4.5) 
r 2 J 3 r 

The equations for 5.F is the same with the one for a global U(l) vortex with the winding number 
1/3. 

First let us solve Eq. (I4.2p . The equation with the right hand side being zero, [A — m| — 
l/(9r 2 )]8F = 0, has a solution q^Ki/^m^r) with being an integration constant. Here Ki/ 3 (r) 
is one of the modified Bessel function K n {r) of the 2nd class, whish solves 

A - m 2 - — KJmr) = 0. (4.6) 
r 2 J 

It is known that K n {r) with < n < 1 is well approximated by Ki/ 2 (r) = y / ^ e_r at r 3> 1. So 
we get 



^^^ e "^H"^ +o ((^ n - (47) 

Since K n {r) is much smaller than 1/r 2 for r ^> 1, the first term can be neglected as in the case 
of the U(l) global (superfluid) vortex. We will not discuss the coefficient below. 



8 If wc consider local vortex by gauging [/(1)b, all the linearized equations at asymptotic region become the 
form of (A — m 2 x )5X = [26]. Here symbol X stands for all the fields, namely Abelian, non-Abelian gauge fields 
and all the scalar fields. 
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Let us next solve Eqs. (14.31) and ( 14. 5ft . Clearly, the solution cannot have a tail decreasing 
with a polynomial. They have exponentially small tails similarly to the ANO vortex in the 
superconductor. The equations with the right hand sides of Eqs. (14.31) and (14.51) being zero can 
be solved by 



SG ~ q x K 1/z {m x r) ~ <?x W ^^e- m * r , Sh ~ q G m G rK x {m G r) ~ q G ^—— e - m ^, (4.8) 

where q x and gc are integration constants, which are determined numerically in the next section. 
Let us take into account the corrections from the right hand sides of Eqs. (14 .3D and ( 14.51) . When 
the inequality rriQ ^ m x holds, we can ignore Sh ~ e~ mGr in Eq. (14.31) . So SG ~ e~ mxr in Eq. (14.81) 
is a good approximation. However, we cannot discard SG in Eq. (14.41) since SG ~ e~ mxT 3> e~ mGT . 
We thus find the following approximations 



7T 



SG = q x < j — 

l x m t% mr for m ^^x(7i^l). (4.9) 

Sh = — ^ „ G „ , /- e~ m * r 



3 — m x Y 2m x r 



On the other hand, when the other inequality m x ^ mc holds, we should reconsider approxima- 
tion for SG. In the same way we find the approximations 



2 <7g \Txm G r _ mar 



3 (™L m c) r2 V 2 6 for m ^ mG ( 7l <i). (4.10) 

Sh~q G J^—e 

These are not good approximation for m x w mg. For the regions m x ^ mg, we find from 
Eqs. (14.91) and (14.101) that both the massive traceless scalars x an d the massive gauge fields 
(gluons) have the same asymptotic behavior with the exponential tails e~ mr with the common 
mass m = min{m x ,mc}. These asymptotic behaviors are quite different from the ANO vortex 
in usual superconductors^ 

Now we can roughly estimate the width of the color magnetic flux tube given in Eq. (13. 6p . 
It can be read from the tail Sh of the gauge fields. Namely it can be estimated by the Compton 



9 Remember that the exponential tails of the ANO vortices are like e~ m " r for the massive scalar field H 
with mass mn and e~ m " r for the massive (7(1) gauge boson with mass m e . The layer structures are exchanged 
for m e > rriH (type II) and m e < mn (type I). Then their interactions depend on the ratio m e /mH giving 
classification of type I and II. However for the strong type II region (m e > 2mjj), the tail of the gauge field 
becomes e~ 2mHr [27]. The width of the gauge field cannot become smaller than the half of that 1/tojj of the 
scalar field. 
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wave lengths m -1 with m = min{m x , mc}. It is the Compton wave length of massive gluons 
(penetration depth) as expected when m x > wiq (7! < 1), while it is no longer the case when 
■me > m x (71 > 1): contrary to the naive expectation the tail of gauge field decays at Compton 
wave length of massive traceless scalar field x (coherence length) l 10 l In Sec. |5] these estimations 
are confirmed by the numerical calculations. 

The U(1)b NG boson makes the power-like tails in the traceless part of $ which is color- 
flavor singlet as seen in Eq. (14. 7p . So the power-like tail is dominated at the large distance from 
the core of the non-Abelian semi-superfluid vortex. It yields the logarithmic divergence in the 
tension of the vortex and leads to the long range repulsive forces between two vortices As 
we approach the vortex core from the large distance, we see the exponential tails behaving as 
e~ m< * r , e~ mxr , e~ mar . We will first come across the tail of the fields with the lightest mass. 

When ?ri0 is the lightest mass, the semi-superfluid vortex behaves as the usual superfluid 
vortex (with 1/3 U(1)b winding), since the massive singlet scalar field <ft exists also in a usual 
£7(1) superfluid vortex. In general, the exchange of the massive scalar field leads to an attractive 
interaction (as in usual superconductors), and so we expect that the short range interaction 
is attractive. In this case non-Abelian properties are somewhat hidden inside the core of the 
superfluid vortex. 

More interesting is the case that m x (or rac) is the smallest, since the asymptotic tail which 
we first come across is those in Eq. (14.91) (or Eq. (14.101) ). There are two remarkable points, i) 
The asymptotic behaviors are quite different from the ANO vortex in usual superconductors. 
The interactions of ANO vortices depend on the ratio of the masses of the massive gauge boson 
and the massive scalar field, giving classification of type I and II as denoted in footnote [9j In 
Eqs. (14. 9 p and (I4.10p we find that both the scalar and vector fields have the asymptotic tails of the 
same order e~ mr with m = min{m x ,mG}, as already pointed out. So we may have to consider 
two contributions equivalently. ii) The interactions may depend on the internal orientation 
moduli CP 2 . In fact, in the case of the local U(N) vortices for which U(1)b is gauged, the 
exchanging of the particles which are not singlet but in the adjoint representation of the color- 
flavor symmetry leads to interactions depending on the internal orientations [28] . In order to get 
better understanding about the short range interactions between non-Abelian semi-superfluid 



10 Similar phenomenon is known for usual superconductors as denoted in the footnote [9] However the threshold 
is m e = 2m h in that case. 
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vortices, we need more qualitative and quantitative studies which is beyond the scope of this 
paper. However see further discussion in Sec. [6j 



5 Numerics: Solutions and Width of Color Magnetic Flux 

In this section we provide numerical solutions. We use the relaxation method with appropriate 
boundary conditions, see Fig. [TJ We see that the interiors of the vortices depend on the ratios 

u, 1, 1) 



1.0 - 

0.8 



(1,5, 1} 



2 4 6 8 10 



{5, 1, 1) 



(1, 1,5} 



{5, 5, 1} 



4 6 
{1,5,5} 



{5, 1, 5} 



02468 10 02468 10 02468 10 

Fig. 1: The numerical solutions for several {big, m <f>: m x\- [M r )> /( r )i 9i r )\ = [R (solid line),G (short dash line),B 
(long dash line)]. 

7i and 72 . 

Some numerical solutions were previously found in [10] with the approximation of g(r) = 1 
in Eq. (13.21) . Although this approximation makes numerical calculation easier, the authors of 
[TO] did not evaluate the validity of the approximation. Note that g(r) = 1 is never true unlike 
the case of non-Abelian global vortex [25J and the non-Abelian local vortex [2B] , in whose cases 

11 



the unwinding fields can be constant for particular values of coupling constants. Moreover, 
G = f — g plays a role of the order parameter for the breaking of the SU(3)c+f, so we need 
the profile function g(r). As we can see from the boundary condition (13.101) . g(0) is not fixed 
and this is a typical parameter to characterize the solutions. We find < g(0) < a/3/2 where 
g(0) — > as m x — > oo and g(0) — > a/3/2 as — > oo. This can be understood by considering 
m x — > oo in Eq. fl 2 . 1 8 [) which forces f 2 + 2g 2 = 3 while — > oo imposes f = g- Qualitatively 
speaking, we find (7(0) < 1 when < m G , and 1 < g(0) for m x < m^. In the case of = m x , 
the profile function g(r) is almost 1 everywhere 11 !. These properties can be seen in Fig. [TJ 

From the solutions in Fig. [TJ we determine the integration constants q x and q G in the asymp- 
totic solutions in Eqs. (14.91) and (14.101) as summarized in Table [TJ 



m G (m x = = 1) 


2 


3 


4 


5 


Qx (SG) 


1.9 (±0.1) 


1.60 (±0.03) 


1.46 (±0.03) 


1.38 (±0.02) 


q x (5h) 


2.0 (±0.1) 


1.62 (±0.01) 


1.47 (±0.02) 


1.38 (±0.01) 


m x (m G = = 1) 


2 


3 


4 


5 


q G (5G) 


1.9 (±0.1) 


1.65 (±0.1) 


1.56 (±0.05) 


1.54 (±0.05) 


q G (5h) 


1.68 (±0.03) 


1.53 (±0.02) 


1.48 (±0.02) 


1.47 (±0.02) 



Table 1: Several numerical values of q x and qa- The case of ma = 2,3,4,5 > 



1 [in Eq. JO])] and the 



case of m x = 2,3,4,5 > mG = 1 [hi Eq. (|4.10[) ] are written in the upper and lower boxes, respectively. The 
parentheses beside q x and qc imply that we evaluate these values by using the the asymptotic tails 8G or Sh in 
Eq. (|4.9p or (|4.10[) . For q x we find a good agreement between q x (SG) and q x {8h). However for qc it is not so 
good for smaller m x . 



Let us next focus on the color magnetic flux by looking at h(r). The function h(r) behaves 

almost similar among the cases of {m G , m^, m x } = {1, 1, 1}, {1, 5, 1}, {1, 1, 5}, {1, 5, 5} in Fig. [1] 

because of the common choice of m G = 1. When the gluon mass is larger, m G = 5, h(r) becomes 

more sharp as can been seen in the cases of {5, 1, 1} and {5, 5, 1} in Fig. [TJ The flux becomes 

far more sharp when we take m x ~ m G as can been seen in the case of {5,1,5} in Fig. [TJ 

Qualitatively these behaviors match the result from the asymptotic tails in the previous section. 

11 The approximation g = 1 used in Ref.|10j may be justified in the region ~ m x . However, since the 
relation Ai = A2 holds at weak coupling regime |10j the masses are related by = 2m x . 
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In order to get a more quantitative width, let us define a width of the color magnetic flux bjl 



The numerical values of this quantity are shown in Fig. [2] for various values of ma and m x . The 
left-upper panel of Fig. [2] shows that the width rapidly decreases for relatively small mc(< 4) 
and then approaches to a constant value for larger mc(> 4). Contrary to this the Compton 
wave length ttiq 1 of the massive gluon goes down to zero as rriQ increases. These behaviors are 
significantly different. Furthermore the width approaches to smaller constant at large mc for 
larger m x as can be seen in the upper-right panel of Fig. [2j We thus find that the size of the color 
magnetic flux cannot become smaller than some value determined by the mass m x of the traceless 
scalar fields for the larger gluon mass ma- This is consistent with the asymptotic behavior found 



On the other hand, the size of the color magnetic flux is well approximated by the Compton 
wave length m^ 1 when vhq changes with keeping m x = viiq as shown in the lower panel in Fig. [21 
The result in this parameter region is complement to the last section, where it was difficult to 
study analytically. 

6 Conclusion and Discussion 

In this paper we have studied the topologically stable vortex in the color superconductor. By 
using the relaxation method with the appropriate boundary conditions, we have presented numer- 
ical solutions for diverse choices of the coupling constants (Fig.[T]). We have found the asymptotic 
tails of the profile functions in Eqs. (14.91) and (I4.10p and have determined their coefficients q x 
and qc in Table [TJ Furthermore we have proposed the width (r) of the color flux tube as in 
Eq. (15 .ip and have calculated it numerically. Contrary to the naive expectation, the width of the 
color-magnetic flux tube does not behave as the Compton wave length (the penetration depth) 
ttiq 1 of the massive gluons for the gluon mass larger than the mass m x of the traceless scalar 
X as shown in Fig. [2j 

12 The definition of this width has been given first in study of a local U(N) vortex for which the phase 
U(1) B is gauged. 




(5.1) 



in Eq. (Q]| . 
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Fig. 2: The width (r) of the color flux tube in Eq. (|5.ip v.s the Compton wave length ttiq 1 . In the upper-left panel 
we plot the width (r) as dots in the cases of m x = 1 and various mo (with = 1). The Compton wave length 
2.2/mc is drawn by the broken curve in which the coefficient is determined by trying to fit it to the numerical 
points. (However they fit well nowhere). The unbroken curve is 2.78 x 7tig~ - 64 made by MATHMATICA which 
fits very well at wiq < 1. In the upper-right panel the width (r) in the cases m x — 1,2,3 (disk, box, diamond) are 
plotted (with = 1). In the lower panel we plot the width (r) in the cases of ttiq — m x (with = 1). The 
broken curve denotes the Compton wave length 2.4/mc which fits to the data better than the case of upper- left 
panel. The unbroken curve denotes 2.72 x m^ ' 82 which fits well to the data everywhere. 
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Here let us discuss the application of our result to the analysis of the interaction between 
two semi-superffuid vortices. As shown in the long range interaction is mediated by the 
massless NG mode of the broken U(1)b- The universal repulsion between two vortices at distance 
R(^$> Tn^ x G ) was found [U] to be 

47T 



F\ ons (R) 



NcR' 



iV c = 3. 



(6.1] 



This force comes from the overlapping of the long power tails of the two vortices. When the 
vortices are placed at a relatively short distance R > m ^ X G^ we ma Y find the exponential tails 
e~ m< * r , e~ mxr , e~ mar . The overlapping of these tails yields the inter-vortex potential. For instance, 
when the tail of x is overlapped for m x < ma, it can be written as 



VJR) = C x 2irqiK {m x R). 



(6.2) 



Here R is the relative distance and q x is the parameter in Eq. (14.91) . The function C x is unknown 
and should depend on the relative orientation of the two vortices in the internal spaceo Thus 
the corresponding inter-vortex force is given by 



9Vx 
dR 



C x 2nq x m x Ki ijn x R) . 



(6.3) 



Let us compare these forces -Fiong an d F x by extrapolating them near the vortex core. In the 
following we assume C x = 0(1) for simplicity. As an example we choose {m x , m^, ma) = (1,1,3), 
see Fig. [3J The width (r) can be read as r = 1.58 from Fig. The approximation G = 



{3, 1, 1} 





Fig. 3: The profile functions for (m x ,m < f > ,mQ) = (1,1,3) in the left-most panel [h(r), f(r), g(r)] =[R (solid 
line),G (short dash line),B (long dash line)]. The middle panel shows the ratio g / (q x Ko(m x r)) ~ 1. A rough 
estimation of the long range and short range interactions are shown in the right-most panel. 



q x Ko(m x r) , initially obtained for the large distance, is still valid up to r > 2 as can be seen in 
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For usual superconductors, C x 



T and this gives the attractive force. 
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the middle panel of Fig. [3j Therefore we can safely use Eq. ( 16. 3ft in this region. As one can see in 
the right panel of Fig. [3], the two kinds of the inter- vortex forces F\ ong and F x extrapolated to the 
near region of the core are comparable around R = 2. We conclude that the inter- vortex force 
mediated by the exchange of the massive particles like the massive gluons becomes comparable 
with the inter-vortex force mediated by the exchange of the massless U(1)b NG boson. More 
precise discussion remains as a future problem. The detailed study of it is necessary to understand 
the lattice structure of non-Abelian semi-superfluid vortices when the lattice spacing is of the 
order of the penetration depth or the coherence length. 

If the short range force is attractive the vortex lattice collapses for high vortex density. 
There may appear composite vortices. The internal structure of a composite vortex may be 
ample which was studied in the case of local U(N) vortices [29J. Other interesting aspects are 
dynamical collision (for instance the reconnection) of non-Abelian vortex strings [30] arid a gas 
of non-Abelian vortices at finite temperature [31] (lower than transition temperature to quark 
gluon plasma phase), both of which were also studied in the case of local U(N) vortices. 

Let us make brief comments on possible application to the neutron star physics. In a response 
to the rapid rotation of a neutron star, superfluid vortices are created along the rotating axis in 
the outer core (the hadronic phase) where neutrons exhibit superfluidity. It was pointed out in 
the seminal paper by Anderson and Itoh [32] that those vortices may explain the pulsar glitch 
phenomenon. If the CFL phase is realized in the core of a neutron star, those superfluid vortices 
are somehow connected to U(1)b superfluid vortices in the CFL phase. Each U(1)b superfluid 
vortex must be divided into three semi-superfluid vortices (color flux tubes). Therefore color 
flux tubes are necessary ingredients in the study of neutron stars if the CFL phase is realized in 
the core. On the other hand, pulsars are accompanied with the large amount of magnetic fields 
of C/(1)em- Protons exhibit superconductivity in the hadronic phase, and so there should exist 
magnetic flux tubes along the axis of the magnetic fields. If those magnetic fields penetrate into 
the CFL phase realized in the core, there is no topological reason for those fluxes to be squeezed, 
because the C/(1)em symmetry (with mixed with one color component) is not broken there. It 
was, however, pointed out that fluxes are squeezed into flux tubes there [13] where the CFL phase 
is modified by large magnetic fields [33]. In general, the axes of the rotation and the magnetic 
fields of the pulsar are different with some angle, the interaction between the (electro-)magnetic 
fluxes and the color magnetic fluxes studied in this paper should be important in the study of 
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the neutron stars. 
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A General Diagonal Solutions 

In this appendix we study more general diagonal vortex solutions. They can be obtained as 
follows: 

B) = «(^-^ nTM) ) fm To _ ^1 {V3 t s + „ 8 r s) ) , (A.i) 

A,(r,«) = -^[\-h(T)\^{^% + ^%), (A.2) 

where we have defined the Cartan subalgebra {T ,T 3 ,T 8 } of U(3) by 

T = -J=diag(l,l,l), T 3 = -^diag(0,l,-l), T 8 = J=diag(-2, 1, 1). (A.3) 

The parameters and v% are determined by the requirement of the single valuedness of the 
fields, eKT^-v^^sTs)) = L Three solutiong are found . ^ jI/g j = (o,l),(±v^/2,-l/2). 

For instance in the case of (z/ 3 , z/ 8 ) = (0, 1) we find the solutions for Eq. (13.11) with 

, F + 2G F-G / . .\ 

/ = -3~. 9 = —. (A.4) 

This solution with (z/ 3 , v 8 ) = (±y3/2, —1/2) constitute the weight vectors of 3. 
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